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RESIZING PROCEDURE FOR OPTIMUM DESIGN 


OF STRUCTURES UNDER COMBINED MECHANICAL 
AND THERMAL LOADING 

By Howard M. Adelman and R. Narayanaswaml 
Langley Research Center 


SUMMARY 


This paper describes an algorithm for resizing structures subjected to 
combined thermal and mechanical loading. The algorithm Is applicable to uni- 
axial stress elements (rods) and membrane biaxial stress members. The proce- 
dure, called thermal fully-stressed design (TFSD) Is based on the basic 
difference between mechanical and thermal stresses In their response to 
resizing. Namely, the mechanical stresses are more sensitive to structure 
resizing than are the thermal stresses. As a result, the TFSD technique Is 
found to converge In fewer Iterations than ordinary fully-stressed design (FSD) 
for problems wh^re thermal stresses are comparable to the mechanical stresses. 

The Improved convergence Is demonstrated by example with a study of a 
simplified wing structure, built-up with rods and membranes and subjected to a 
combination of mechanical loads and a three-dimensional temperature distribu- 
tion. Both the FSD and TFSD methods converged to the same final design and 
TFSD required far fewer Iterations to converge than did FSD. 


INTRODUCTION 


Probably the most widely used approach for sizing of flight structures 
under strength and minimum gage constraints Is fully-stressed design (FSD). 

In this method the structural sizes are Iterated with the step size depending 
on the ratio of the total stress to the allowable stress (refs. 1-3). The FSD 
procedure Is traditionally used to obtain, at a reasonable computational cost, 
designs which If not at a minimum weight are at least acceptably close to the 
minimum weight (ref. 2). 

Almost all of the experience with FSD has been with structures primarily 
under mechanical loading as opposed to thermal loading. The temptation In 
Including thermal loads In FSD Is to simply continue to use the total stresses 
In computing the Iteration step size. This approach seems satisfactory when 
mechanical stresses dominate the thermal stresses (ref. 4). Convergence may 
be slo\.’, however, when thermal stresses are comparable to mechanical stresses. 
The slower convergence Is associated with relative Insensitivity of the thermal 
stresses to changes In structural sizing. Procedures are therefore needed 
which take Into account the differing responses of thermal and mechanical 
stresses to changes In structural sizes. 


An Improved variant of FSD was described In reference 4 for uniaxial 
stress members. It was demonstrated for automated sizing of truss-type struc- 
tures. For problems having substantial thermal stress • the new procedure 
(called thermal fully-stressed design or TFSD) was found to converge in far 
fewer Iterations than ordinary FSD. This paper extends the TFSD procedure to 
biaxial stress members using the Von Mlses failure criterion. The TFSD resizing 
procedure for uniaxial stress members Is restated^ the new procedure for biaxial 
stress members Is developed, and results are given from an application of the 
procedure to size a simplified wing structure. 
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LIST OF SYMBOLS 
cross sectional area of a bar 
thermal-mechanical stress coupling term (eqn 8) 

Young's modulus 

number of iterations required for FSD to converge to within 5 percent 
of final mass 

number of Iterations required for TFSD to converge to within 5 per- 
cent of final mass 

” ^i+l^*^i 


thickness of a membrane 

■ p + - a a + 1 , 

lx y XV xyj 


Von Mlses stress measure 


coefficient of linear thermal expansion 
weight density 
Poisson's ratio 
stress component 


Subscripts 

a allowable 

i Iteration number 

M mechanical 

T thermal 

X, y orthogonal coordinate directions in plane of membrane element 
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THE TFSD ALGORITHM FOR UNIAXIAL STRESS MEMBERS 


The TFSD resizing algorithm for uniaxial stress members (rods) Is given In 
reference 4 as 


A 


1+1 



( 1 ) 


In equation (1) Is the stress due to mechanical loads acting alone, 

is the stress due to thermal loads acting alone and a is either the 

a, M 

tensile or compressive allowable stress, depending on the sign of a^. The 
algorithm of equation (1) drives each element toward the condition 
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„ Cl ^ 
a,M - T 


- 1 


( 2 ) 


Thus, the mechanical stress Is driven toward an effective allowable stress 
given by the algebraic difference between the material allowable stress and 
thermal stress. For later reference, the usual FSD algorithm Is given by 
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(3) 


where in equation (3), 0 Is either the tensile or compressive allowable 

d 

stress, depending on the sign of the total stress as given by the numerator. 
The algorithm In equation (3) drives each element toward the condition 
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DEVELOPMENT OF TFSD RESIZING ALGORITHM 
FOR BIAXIAL STRESS MEMBERS 


The Von Mlses failure criterion for Isotropic biaxial stress members is 


3 
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- 0 


( 4 ) 


where 0^ and are direct stresses along orthogonal coordinate directions 

In the plane of the element, 0^ Is the shear stress on that plane and 

Is the allowable stress. In preparation for generalizing the resizing formula 
from uniaxial to biaxial stress members, rearrange equation (1) as follows: 


'Ml 
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(5) 


By analogy, the corresponding statement for blaxlally stressed members Is 
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( 6 ) 


where ■ ^1+l^^i element thickness. Expansion of equation (6) 

using equation (4) gives 


Rl - “a) ■^1 " "I'^i * <1 • “ 


(7) 


where 


9 7 2 2 

V ^o +0 -a^o^+3o _ 

'^1 xT yT xT yT xyT 
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( 8 ) 



Solution of equation (7) by the quadratic formula gives the resizing algorithm 



The choice of sign In front of the radical Is dictated by the requirement that 
the bracketed quantity be positive since only positive thicknesses are physi- 
cally meaningful. When the positive sign must be chosen. When 

V > a and b > 0 the algorithm Is Inapplicable since no acceptable design 

X d 

exists without thermal stress reduction by non-structural means. When 

V > a and b < 0, either the positive or negative sign may be used provided 

X A 

the radlcand Is positive. In the later case Intuition would suggest the choice 
of the negative sign since the minimum weight design Is sought. The authors' 
use of the algorithm has been limited to cases wherein V < o and hence the 

X d 

positive sign was used exclusively In equation (9). 

For later reference, the corresponding FSD resizing algorithm used to 
compare results with TFST) Is 


t 


1+1 



( 10 ) 


where V is defined by equation (4). 


APPLICATION TO A BUILT-UP WING 


To Illustrate the application of the TFSD resizing algorithm and to com- 
pare the algorithm with ordinary FSD, calculations were carried out using a 
computer program incorporating both the TFSD and the FSD procedures. Finite 
element methods using standard rod elements and the "TRIM 6" (ref. 5) triangu- 
lar membrane elements were used for the analyses. 

The structure used in the calculations is a simplified low-aspect ratio 
built-up wing structure (shown in fig. 1) which is roughly based on the configu- 
ration studied extensively in reference 6. The ribs and spars are modeled by 
trusses with a total of 85 rod elements. The upper and lower skins are each 
modeled by 5 membrane elements. As a result there are a total of 95 design 
variables in the problem including rod areas and membrane thicknesses and the 
finite element model has 36 grid points. The finite element model is described 
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In Table 1 where grid point coordinates are listed and in Table 2 which gives 
the connection grid points of each rod element. The loads on the wing consist 
of concentrated forces representing pressure loads of 13.8 KPa (2 psi) in the 
positive Z - direction at points on the lower surface in addition to elevon 
loads and a three-dimensional temperature distribution in which temperatures 
range up to about 600 K. The loads are tabulated in Table 3 and the tempera- 
ture distribution is shown in figure 2. 

Identical results were obtained for the final design by the TFSD procedure 
(equations (1) and (9)) and by the FSD procedure (equations (3) and (10)). 

The final design for the rod areas and membrane thicknesses are given in Tables 
4 and 5. respectively. As an aid in Interpreting the final design, a pictorial 
representation of the distribution of membrane thickness is shown in figure 3. 
The largest thicknesses on both the upper and lower surfaces are in the regions 
of the trailing edge of the wing, namely elements 4 and 5 on the upper surface 
and elements 9 and 10 on the lower surface. In these regions, the in-plane 
mechanical loads are maximum due to the presence of the elevon loads, while the 
thermal loads which relieve the mechanical loads are smallest. As a result, 
the lar>;;< St total stresses occur in membranes 4, 5, 9, and 10. 

The Important observation from this calculation is made by comparing the 
number of des^ n iterations required to converge from an arbitrary trial design 
in which all ileslgn variables had unit values. The FSD algorithm required 18 
iterations to converge to within 5 percent of the mass of the final design 
whereas the TFSD obtained this degree of convergence in a single iteration. 

This type of performance by the TFSD procedure relative to FSD is consistent 
with previous experience comparing the two methods. In reference 4, TFSD 
required one-fourth as many iterations as FSD to converge for a structure 
modeled by rod elements only. 

To further investigate the convergence of TFSD relative to FSD, three addi- 
tional sets of calculations were performed; this time with lower temperature 
levels but with the same mechanical loads as in Table 3. The TFSD procedure 
required only a single iteration to converge to within 5 percent of the final 
mass for all cases. A summary of the convergence of TFSD relative to FSD is 

shown in Table 6 and figure 4 as a function of T/T . where T . represents 

ret ret 

the highest level of thermal loads and corresponds to the temperatures in 
figure 2 and Table 3. In the first of these cases, the temperatures were input 
as 75 percent of those in Table 3 “ *75). The design for this case is 

given in Tables 7 and 8 and is 6 percent lighter than the design for the refer- 
ence temperatures. The FSD algorithm required 10 iterations to converge to 
within 5 percent of the final mass. In the second case the temperatures were 
30 percenc of those in Table 3. The final design is given in Tables 9 and 10 
and is 9 percent lighter than the design corresponding to the reference temper- 
atures. FSD required 5 iterations to converge within 5 percent of the final 
mass. In the third case, the temperatures were 25 percent of those in Table 3. 
The final design given in Tables 11 and 12 is 10 percent lighter than the design 
for the reference temperatures and FSD required 4 iterations to converge to 
within 5 percent of final mass. 
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Table 6. - Effect of Temperature on Relative Efficiency of 
FSD And TFSD for Sample Problem 


I/Tr„ 

0. 

.25 

.50 

.75 

1.0 

^FSD 

1 

4 

5 

10 

18 

V 

TFSD 

1 

1 

1 

1 

1 


For this example as well aa for the examples of reference A, thermal 
stresses are quite insensitive to structural sizing. The superiority of TFSD 
for these examples is associated with this insensitivity. This will be clear 
if we imagine a case where thermal stressea are completely independent of struc- 
tural size in which case TFSD would obviously be superior. It seems reasonable 
that a broad range of structures will exhibit relative insensitivity of thermal 
stresses to sizing. Consequently, the TFSD procedure should be widely useful 
for structures under combined thermal and mechanical loads. 


CONCLUDING REMARKS 


The paper presents an Improved algorithm for resizing structures subjected 
to combined thermal and mechanical loading. The algorithm originally developed 
for resizing uniaxial stress elements (rods) is herein extended to biaxial 
stress elements such as membranes. The algorithm is based on monitoring the 
mechanical and thermal stresses separately and altering the mechanical stresses 
rather than the total stresses. The thermal stresses enter into the algorithm 
by adding to or subtracting from the allowable stresses. The improved algo- 
rithm, called thermal fully-stressed design (TFSD) takes account of the basic 
difference between mechanical and thermal stresses in their response to 
resizing. Namely, the mechanical stresses are more sensitive to structure 
resizing than are the thermal stresses. As a result, the TFSD technique tends 
to converge in fewer Iterations than does ordinary fully-stressed design (FSD) 
for problems where thermal stresses are comparable to the mechanical stresses. 

This behavior is demonstrated in the paper by calculations on a simplified 
built-up wing under a combination of mechanical loads and a three-dimensional 
temperature distribution. Both the FSD and the TFSD methods converged to the 
same final design, but TFSD converged much more rapidly than did FSD. Further, 
the relative efficiency of TFSD relative to FSD increased according to the level 
of the applied thermal loading. 
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TABLE I. - CCX)RDINATES OF GRIDPOINTS 
OF WING FINITE ELEMENT MODEL 


GRID 

X 

Y 

Z 

POINT 

cm 

in 

cm 

In 

cm 

in 

1 

0 


0 

0 



12 

2 

G 


0 

0 

0 


0 

3 

190.5 

75 

0 

0 


12 

4 



75 

- 50.8 

-20 


12 

5 


1 

75 

0 

0 


0 

6 



75 

- 50.8 

-20 


0 

•/ 

381.0 

150 

0 

0 


12 

3 



150 

- 50.8 

-20 

30.48 

12 

9 



150 

-101.6 

-40 

30.48 

12 

10 



150 

0 

0 

0 

0 

11 


1 

150 

- 50.8 

-20 

0 

0 

12 



150 

-101.6 

-40 

0 

0 

13 

571.5 

225 

0 

0 

30.48 

12 

14 



225 

-101.6 

-40 

30.48 

12 

15 



225 

-152.4 

-60 

30.48 

12 

16 



225 

0 

0 

0 

0 

17 



225 

-101.6 

-40 

0 

0 

18 



225 

-152.4 

-60 

0 

0 

19 

762.0 

300 

0 

0 

30.48 

12 I 

20 



300 

-101.6 

-40 

30.48 


21 



300 

-203.2 

-80 

30.48 

12 1 

I 22 



300 

0 

0 

0 

0 1 

1 23 



300 

-101.6 

-40 

0 

0 1 

; 24 



300 

-203.2 

-80 

0 

0 ; 

: 25 

825.5 

325 

0 

0 

30.48 

12 ! 

i 26 



325 

-101.6 

-40 

30.48 

12 : 

1 27 



325 

-203.2 

-80 

30.48 

12 : 

! 28 



325 

0 

0 

0 

0 ! 

I 29 

1 


325 

•101.6 

-40 

0 

0 ; 

30 



325 

-203.2 

-80 

0 

0 ! 

31 

889.0 

350 

0 

0 

30.48 

12 ! 

32 



350 

-101.6 

-40 

30.48 

12 1 

I 33 



350 

-203.2 

-80 

30.48 

12 

: 34 



350 

0 

0 

0 

0 

35 



350 

-101.6 

-40 

0 

0 : 

36 



350 

-203.6 

-80 

0 
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TABLE 2 


CONNECTIVITY TABLE FOR 























TABLE 4. - FINAL DES?GN OF TRUSS ELEMENTS OBTAINED BY TFSD AND FSD 
Tenq>eratures and Loads Given In Tablo 3. 
















TABLE 5, - FINAL DESIGN OF SKIN (MEMBRANE) 
ELEMENTS BY FSD AND TFSD METHODS 

Loads and Temperatures Given in Table 3. 


Element 

Number 

r 

Thickness 

cm 

in 

1 

.002540 

.0010 

2 

.006812 

.002683 

3 ' 

1 

.002540 

.0010 

^ ; 

1 

.03383 

.01332 

= 1 

.03180 

.01252 

^ i 

.002540 

.0010 

^ i 

.02195 

.008640 

« 1 

1 

.005022 

.001977 i 

1 

9 

1 

.02652 

! 

.01044 [ 

! 10 i 

1 

1 

1 

.03160 

.01244 


Final Mass 30.78 kg (67.80 Ibm) 


TABLE 7 


FINAL DESIGN OF TRUSS ELEMENTS OBTAINED BY TFSD AND FSD METHODS 


Temperatures are 75 Percent of Those in Table 3. 


Bar 

— 

A 

rea 

Bar 

Area 

Bar 

Area 

2 

I cm 

— 

in^ 

2 

cm 

in^ 

2 

cm 

in^ 

1 

.00( 

54 

.001 


30 

) 

• 

0064 

.001 

59 

.0064 

.001 

2 





1 




1 

60 




1 

3 





t 32 





61 



1 


4 

' 




1 33 





62 

.3034 

.04703 

5 

.04543 

.007042 

i 34 




I 

1 

63 

.0064 

.OC 

1 

6 

.0064 

.001 

35 





64 





7 

.0064 

.001 

36 




i 

65 





8 

.06813 

.01056 

37 




1 

1 

66 





9 

.0064 

.001 

: 38 




i 

I 

67 





10 

.0064 

.001 

1 

1 

1 

39 




, 

68 

.6065 

.09401 

11 

.09845 

.01526 

' 40 

1.315 

.2038 1 

69 

.0064 

.001 

12 

.0064 

.001 

1 

I 41 

.0064 

.001 

70 





13 ! 

.0064 

.001 

1 

1 

1 

42 




1 

1 

1 

71 





14 

.04543 

.007042 

43 




1 

\ 

72 


1 

I 



'' i 

.0064 

.001 

1 

1 

1 

44 





73 





16 1 

.0064 

.001 

1 

1 

45 





74 





17 

.01514 

.002347 i 

46 





75 


t 


1 

18 

.0064 

.001 

1 

47 





76 


j 



19 

.132 

4 : 

.02053 ' 

48 





i 

77 


1 



20 

.08832 

.01369 : 

49 





78 j 


1 

1 



21 

.353 

0 i 

i 

.05475 ' 

\ 

50 





t 

79 1 

1 

.1033 ' 

.01601 

22 

.1324 

.02053 1 

31 

.8593 

.1332 

80 1 

.0064 , 

.001 

23 

.9484 

.1470 j 

52 ; 

.0064 

.001 

81 


1 

1 



24 

. 342 

6 1 
1 

.05310 i 

53 , 





82 


1 

t 

1 



1 25 

1 

.974 


.1511 i 

54 i 


1 



83 





26 

.1581 

.02451 1 

55 ! 





84 ‘ 





27 

.4743 

.07352 ' 

56 





85 

.5634 1 

.08733 

28 

.052 

71 

.008169 1 

57 i 




1 






; 

.2108 j 

.03268 j 

58 

— • ^ ^ ^ J 




1 





J 


14 
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TABLE 8. - FINAL DESIGN OF SKIN (MEMBRANE) 
ELEMENTS BY TFSD AND FSD METHODS 

Temperatures are 75 Percent of Those In Table 3. 


Element 

Number 

Thickness 

cm 

In 

1 

.00254 

.001 

2 

.007313 

.002879 

3 

.00254 

.001 

4 

.03165 

.01246 

5 

.03096 

.01219 

6 

.00254 

.001 

7 

.01639 

.006452 

8 

.003325 

.001309 

9 

.02670 

.01051 

10 

.03193 

.01257 


Final Mass 29.10 kg (64.01 Ibm) 



TABLE 9. - FINAL DESIGN OF TRUSS ELEMENTS OBTAINED BY TFSD AND FSD METHODS 


Temperatures Jne-half Values in Table 3. 





















TABLE 10. - FINAL DESIGN OF SkTIN (MEMBRANE) 
ELEMENTS BY TFSD AND FSD METHODS 
Temperatures One-Half Values in Table 3. 


Element 

Number 

Thickness 

cm 

in 

1 

.002540 

.001 

2 

.007714 

.003137 

3 

.002540 

.001 

4 

.02964 

.01167 , 

! 

5 

.03015 

> 

.01187 

• 

6 

.002540 

.001 

7 

0 

.01218 

.004796 

8 

.002540 

.001 

9 

.02814 

.01108 

10 

.03210 

L- - 

.01264 


Final Mass 28.09 kg (61.80 Ibm) 


TABLE 11. - FINAL DESIGN OF TRUSS ELEMENTS OBTAINED BY TFSD AND FSD METHODS 


Temperatures are One- fourth Values in Table 3. 


Bar 

Area 

— 

Area 

Bar 

Area 

2 

cm 

in^ 

Bar 

2 

cm 

In^ 

2 

cm 

in^ 

1 

.00( 

S4 

.001 

30 

1 

• 

0064 

.0 


59 

.0 

064 

.00 

1 

2 





' 31 




1 

60 





3 





' 32 

1 




( 

61 





A 





1 33 





62 

.3494 

.05415 

5 

.04543 

.007042 

34 





63 

.0064 

.00 

1 

5 

.0064 

.001 

35 





64 





7 

.0064 

.001 

1 36 





65 





8 

.06813 

.01056 

' 37 

1 





66 





9 

.0064 

.001 

j 38 





67 





10 

.0064 

.001 

39 

1 





68 

.7445 

.1154 

11 

.09845 

.01526 
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TABLE 12. - FINAL DESIGN OF SKIN (MEMBRANE) 
ELEMENTS BY TFSD AND FSD METHODS 

Temperatures One— Fourth Values in Table 3. 


Element 

Thickness 

Number 

cm 

in 

1 

.002540 

.001 

2 

.008304 

.003466 

3 

.002540 

.001 

4 

.02776 

.01093 

5 

.02941 

.01158 

6 

.002540 

.001 

7 

.009030 

.003555 

8 

.002540 

.001 

9 

.03007 

.01184 

10 

.03216 

.01266 

1 




Final Mass 27.57 kg (60.65 Ibm) 




(a) Grid point locations and Truss elements 



Upper skin 

(b) Membrane elements 


Lower skin 


Figure 1.- Finite element model of built-up wing. E = 200 GPa (3 x 10^ psi) 
n = 8303 kg/m^ (0.3 Ibm/in.^), a = 12.0 x 10"®/K (6.5 X 10"®/°F), 
- 0.97 GPa (142 psi), i/ = 0.3. 



433 



(a) Upper surface 


Fij^ure 2.- Temperature distribution on wing. (Values are in Kelvin) 




Fif^ure 2.- Concluded. (Values are in Kelvin) 



FSD 


FSD 


Figure 4.- Effect of thermal load level on convei^encc of FSD as comoared to TP«Jn 


